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1. Introduction 

I review the status of light hadron spectroscopy from lattice QCD. I start by reviewing lattice 
calculations of the masses of flavour singlet pseudoscalar mesons. In the first part I will focus on 
lattice results for the J^'-- = 0++ scalar mesons, because these mesons are candidates for having 
glueball or tetraquark degrees of freedom and are still not well understood. In the second half of 
the talk I will compare the results for basic quantities such as the pion decay constant and mass of 
the light vector meson between different lattice formalisms. This is an important part of validating 
lattice calculations. 

There are various omissions in this review. I don't include any results for baryons. Of par- 
ticular note is work of the LHPC collaboration [|l|] who are using a highly developed variational 
technique to try to fully map out the low lying baryon spectrum. Also I don't discuss any devel- 
opments in the spectroscopy of mesons that include heavy quarks, although there have been many 
new states that have been discovered such as the 0^(2317), X(3872), and F(4260). See for a 
review of the experiments and the results from model calculations, and ^ for a review of lattice 
results. 

2. The singlet pseudoscalar mesons from lattice qcd 

The large mass of the rj' meson is thought to be caused by the QCD vacuum structure and the 
axial anomaly. The rj and rj' mesons are decay products of flavour non-singlet 0++ mesons, so are 
a natural starting point for a review with a focus on 0++ mesons. 

In unquenched QCD with «/ =2 sea quarks there is only one singlet pseudoscalar meson that 
I will use the T]2 notation for. The mass of the t]2 is expected to be around 800 MeV. The lattice 
calculations of the flavour singlet pseudoscalar mesons involve the computation of disconnected 
diagrams, that are more noisy and compute intensive, than connected correlators. With the avail- 
able computing power it is possible to compute the relevant disconnected correlators. However, 
the intrinsic noisiness of the correlators means that many more configurations are required than for 
a standard lattice QCD calculation of flavour non-singlet quantities, and this makes them expen- 
sive [Q, ^. In fact the sub-groups who work on flavour singlet quantities inside a collaboration are 
always the people who ask for much longer simulation runs. 

In figure [l] I show a plot of recent results from the ETM collaboration for the mass of the T]2 
meson [Q]. There is an summary of the results of older lattice QCD calculations for the mass of the 
T]2 meson in [^. 

Figure |l] shows that the mass of T]2 meson is consistent with a constant behaviour with quark 
mass, although the statistical errors are large. The preliminary estimate for the mass of the 772, from 
ETMC [|], is « .88 GeV (rom(T]2) = 2). 

Another complication of 2 + 1 lattice QCD calculations is rj-rj' mixing. In lattice QCD calcu- 
lations with 2-1-1 flavours of sea quarks there will be mixing between light and strange interpolating 
operators. The experience with the last generation of «/=2 lattice QCD calculations with pion 
mass at the 600 MeV level found that the differences between quenched and unquenched QCD 
were mostly small [Q]. This suggests that the difference between «/=2 and rif =2+1 will be small. 
However this will not be the case for the flavour singlet pseudoscalar mesons where the ground 
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Figure 1: Summary of results for the mass of 772 from calculations with light pions and a < 0.1 fm. The 
results are from ETMC [1, CP-PACS [pi], and UKQCD Pn. 



state will be the mass of the rj (mass 548 MeV) with ?iy = 2+1 flavours, but the mass of the 772 
(mass around 800 MeV) for lattice calculations with «/ = 2 flavours. This will indirectly effect 
particles that decay via the strong interaction with a flavour singlet pseudoscalar meson as a final 
state, because the decay threshold will differ between «/ = 2 and = 2 + 1, by the order of 250 
MeV. 

Both the ujsu + dy^d and ly^s will couple to the r\ meson. The r\' is the first excited state in 
light flavour singlet pseudoscalar channels. A modern approach to 77, 77' mixing is reviewed by 
Feldmann [Q]. 

The CP-PACS/JLQCD collaboration used a variational technique to study the 77 and 77' me- 



son [10]. The basis states in equation 2.1 were used to form a variational smearing matrix 2.2 



ri„ = {uY5U + dY5d)/V2, ri, = {sY5S 



(2.1) 



^(0 = I y^^X^ . 3 1 > (2-2) 



The preliminary results presented at lattice 2006 [[To|]. from a lattice calculation with a ~ 0.12 
fm, mv/mps 0.61 to 0.78 were = 0.55(2) GeV, and = 0.87(5) GeV (compare to with ex- 
periment, rriri = 0.548 GeV, and m,,/ = 0.958 GeV There has also been a recent attempt to compute 
the masses of the 77, 77' mesons using improved staggered fermions [^, that is also regarded as an 



important theoretical test of the staggered formalism [ 11 1 



3. The light scalar mesons from lattice QCD 

The interpretation of many 0++ mesons in terms of quarks and glue degrees of freedom is 
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Meson 


I 


MMeV 


TMeV 


main decay 


comment 


ao(980) 


1 


985 


50 - 100 


rjn 


tetraquark, molecule, qq 


ao(1450) 


1 


1474 


265 


r\n 


qq 


K 


1/2 


660 


560 


Kn 


not in PDG (yet) 


^o(1430) 


1/2 


1414 


290 


Kk 


qs 



Table 1: Light 0++ and 0+ flavour non-singlet mesons 



Group 




niKo GeV 


Prelovsek et al. [p^] 


2 


1.6±0.2 


McNeile and Michael [|l^] 


2 


1.1-1.2 


Mathur et al. ^ 





1.41 ±0.12 


SCALAR [0] 





~ 1.7 



Table 2: Lightest strange-light 0+ meson from lattice QCD. 



Still not clear. The 0++ mesons potentially contain glueball, tetraquark, meson molecule or even 
quark-antiquark degrees of freedom. 

There are a number of reasons that lattice calculations of the light scalar mesons are challeng- 
ing. The lattice QCD correlators for scalar mesons are more noisy than for p and n mesons. The 
light scalar mesons decay via S-wave decays, so the state decays when its mass equals the sum of 
the masses of the decay products. This makes it easier to see the effect of the strong decay in lattice 
calculations, than for states that decay via P-wave, such as p. A, or Roper resonance. 



3.1 The flavour non-singlet 0++ and 0+ mesons. 

In table [T] I collect some pertinent experimental properties of the light flavour non-singlet scalar 
mesons. For more detailed information see the reviews [0, 13]. The existence of the K meson is 
controversial (see [ p^ ] for a discussion), I use the masses and widths from [[l^]. 

Although I am going to loop through the lattice results for the different mesons in table [I one 
by one, it is also important to classify the states into SU3 multiplets or a classification based on 
tetraquarks. The key questions we want to answer from lattice QCD are: do we see ao(980) with 
qq operators (where <7 is a generic light quark) and do we see the K meson at all? 

In table ^ I collect results for the mass of the lightest 0+ qs meson from lattice QCD calcu- 
lations. The lattice results in table ^ are consistent with experimental mass of the Kq{\430), but 
mostly miss the controversial K particle. All the lattice calculations used sq interpolating oper- 
ators, so may have missed the K state, if it is mostly a tetraquark state, with no overlap with sq 
interpolating operators. 

I now discuss the mass of the lightest flavour non-singlet 0++ meson from quenched and 
unquenched QCD. In quenched QCD there is a ghost contribution [|^], due to the T] tt contribution, 
to the scalar correlator that needs to be subtracted off the lattice data. This contribution is included 
via chiral perturbation theory jr9|]. I collect together some recent results for the mass of the light 
0++ meson from lattice QCD in table |3[ I only include quenched data where the r] n contribution 
has been corrected, hence there is no data before the paper by Bardeen at al. [ 19] in 2001. 
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Group 


iif 


Mao GeV 


Bardeen at al. ^ 

Burch et al. 
Hart etal. [21] 
Prelovsek et al. [|l|] 
Prelovsek et al. [|l|] 

Mathur et al. [18] 





2P 

2 

2P 



1.34(9) 
~ 1.45 
1.0(2) 
1.58(34) 
1.51(19) 
1.42(13) 



Table 3: A collection of results from lattice QCD for the mass of the lightest non-singlet 0++ meson. The P 
stands for partially quenched. 
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Figure 2: Difference in mass between masses of light bi and ao- 



McNeile and Michael []17p, in an unquenched lattice QCD calculation focused on the mass 
difference (in the hope that systematics cancel), between the 1+^ and the 0^+ mesons. The lat- 
tice calculation used gauge configurations from UKQCD's non-perturbatively improved clover ac- 
tion and configurations from CP-PACS's tadpole improved clover program [22]. The results 
for the mass difference are plotted in figure ^. Figure ^ shows that the mass of the 1+^ was higher 
than the 0++. The final result was mh^ — m„g = 221 (40) MeV, compared to the experimental result 
of 245 MeV. At this conference Lang et al. reported masses for the lightest flavour non-singlet 0++ 
consistent with the mass of the fl;o(980) meson, from an unquenched lattice QCD calculation using 
chirally improved fermions [p3|]. 
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Figure 3: Mass of lightest state in 0++ channel with nrj2 decay threshold. 



The previous lattice QCD calculations were in a regime where the quark masses were large 
enough that the decay oq rjTi was forbidden. Now I discuss the new lattice QCD calculations 
where the decay oq ^ rjn is energetically allowed. 

The MILC collaboration [ 24 ] originally claimed that they had evidence for ao decay to nrj 
from their 2+1 calculations with improved staggered fermions. Other decays are discussed in [^]. 
Later work by the MILC [25] and UKQCD [^] collaborations showed that the lightest state in the 
flavour non-singlet 0++ channel was actually below the nrj threshold, with improved staggered 
fermions. This was puzzhng, because experimentally the ao Tin decay is forbidden by G parity. 

In [27], Prelovsek explained the behaviour of the flavour non-singlet 0++ correlator with im- 
proved staggered fermions using staggered chiral perturbation theory. Bernard, DeTar, Fu, and 
Prelovsek [28] extended the original analysis by Prelovsek, and also applied it to the flavour singlet 
/o meson. This is a successful theoretical test of the rooting of the staggered determinant, but a 
larger study, with more sea quark masses, is required to say something specific about the mass of 
the ao meson. 

The ETM collaboration have preUminary results for the mass of the light 0++ meson from a 
n f=2 unquenched lattice QCD calculation with twisted mass fermions. In figure || I plot the mass 
of the light 0++ meson and the 71 + 172 decay threshold as a function of the square of the pion mass. 
The mass of the 772 was computed by Michael and Urbach [Q]. Figure |3| shows some evidence for 
the mass of the 0++ tracking the 71 + T]2 threshold, or at least for it being an open decay channel. 
There was an error with the original preliminary data from ETMC for the masses of the ao meson 
in figure ^, so the plot has been replaced with the correct one from arXiv:0906.4720. 

Some caution is required in interpreting the results, because we are only just starting to deal 
with mesons with open decay channels in unquenched lattice QCD calculations. As discussed in 
section there is a big difference between the lightest flavour singlet pseudoscalar meson in lattice 
QCD calculations with «/ = 2 and rif = 2 + 1 sea quark flavours. 

To learn how to deal with mesons with open decays on the lattice we need some simple test 
cases to validate the lattice methods. The Zji(1235) meson is good example, because most people 
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Figure 4: Mass of bi state as a function of square of pion mass 



think that it is qq state and its width is not too large at 142 MeV [p9[]. A bad example to study is 
the fli (1260) with an experimental width of 250 to 600 MeV [|^]. The bi meson has the dominant 
decay con, but I will plot the pn decay threshold because the difference between the p and co 
correlators are disconnected and are thought to be small. 

In figure ^ I plot some preliminary results from the ETM collaboration for the mass of the b\ 
meson with the estimate of the con threshold, as a function of the square of the pion mass. There 
was an error with the original preliminary data from ETMC for the masses of the bi meson in 
figure^ so the plot has been replaced with the correct one from arXiv:0906.4720. 

3.2 Determining the structure of a meson 

Some questions such "is this state a molecule?" can be reformulated as: "what is the structure 
of this hadron?". Most lattice QCD calculations just introduce a particular class of interpolating 
operator, such as qq or qqqq and then check to see what the lowest masses are. It would be useful 
if lattice QCD calculations also tried to look at the structure of mesons. This is after all how 
continuum phenomenology is done or attempted. 

In the past there have been attempts to measure wave functions of hadrons of the form \j7{x + 
r)Y{x). This can be done by fixing the gauge, or inserting gauge links between the quark and 
antiquark fields. Because the \j7{x + r)Y{x) operator is restricted to one particular Fock space 
component these objects are not directly accessible to experiment via a form factor measurement. 



but are still a valuable theoretical tool. The CLQCD collaboration [30] have recently argued that 
because they found a node in the lattice wave-function of the first excited 1++ state (a candidate 
for the X(3872)) in charmonium, then the state was a conventional cc meson. 

The one test that is commonly used in lattice QCD calculations, that include tetraquark degrees 
of freedom, is to study the volume dependence of the amplitudes. A strong volume dependence 



indicates that the state is a scattering state and no resonance is formed [31]. 

There are also observables that are sensitive to the structure of the hadron that are also in 
principle accessible to experiment. For example, leptonic decay constants, strong decay widths. 
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Group 


Method 


fao MeV 


Maltman [38] 
Shakin and Wang [39] 
Narison [Q 


sum rule 
model 
sum rule 


298 
433 
320-390 



Table 4: Some results from models for fa^ decay constant. 



and two photon decay widths. 

Pennington [ p2| ] has recently extracted the two photon decay width of the a from experiment 
to be r(a 77) ~ 4 keV. Pennington notes that value of r(a 77) can depend quite sensitively 



on the quark content of the a [32]. Thanks to Dudek and Edwards [B3] we now have a formalism 



to compute two photon widths on the lattice. Dudek and Edwards [ |33[ ] compute r(;^o 77) = 
2.4 zb 1 .0 keV, from a quenched QCD calculation. It would be interesting to do a similar calculation 
for light scalars. 



Narison [34] proposed to use the leptonic decay constant of the non-singlet 0++ mesons to 
determine the structure of the ao meson. The fag decay constant of the light flavour non-singlet 
0++ meson has been computed using unquenched lattice QCD [p^]. 

{0\qq\ao) =MaJaQ (3.1) 



See [ 17] for a further discussion of this decay constant and the connection with the electroweak cur- 



rent. The experimental measurement of the decay constant in equation via T decay is discussed 



by Diehland Hiller [35] 



A molecule of two mesons should have a very small "wave-function" at the origin, hence fa^ 
should be small. The definition of fa^ is similar to that of the pion decay constant. Hence we 
mean "small" relative to 130 MeV. The other measured decay constants of pseudoscalar mesons 
are within a factor of 2 to the pion decay constant [|29|]. The only exception is the decay constant 



of the 71 (1300) that is suppressed | [3q , 37]. A large value for decay constant fa^ does not rule out a 
qqqq multi-quark meson. 

Using gauge configurations from UKQCD and CP-PACS, McNeile and Michael computed 
fao 480 MeV. In table ^ I collect some results for the value of fag from model calculations. A 
lattice QCD calculation of the decay constant of the scalar heavy-light mesons was reported in [pO|]. 

Computing the decay width of a hadron is also very a valuable way of identifying a state on the 



lattice. In [17], it was reported that the experimental hadron coupling for the decays ao(980) KK 
and ao(1450) KK were 0.9 and 0.5 respectively. A lattice calculation [ ]l7| ] found that the lightest 
hadron in the 0++ correlator had a coupling to KK of « 1 , thus providing additional evidence that 
the lightest state was the ao(980). 

3.3 Flavour singlet 0++ mesons 



In pure SU3 gauge theory, Morningstar and Peardon [ |41| ] found 13 glueballs with masses 
under 4 GeV This raises the question as to whether there is any evidence for glueballs in nature, 
or are glueballs just theoretical constructs that are only of interest to test ADS/CFT techniques for 
example. One way we can answer this question on the lattice is to study the effect of sea quark 
dynamics on the glueball masses. I will focus on the 0++ state, as this is where the most work has 
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Meson 


MMeV 


TMeV 


decay 


comment 


/o(600) (a) 


441 


544 


KK 


tetraquark, molecule, qq 


/o(980) 


980 


40 - 100 


71 71 


tetraquark, molecule, qq 


/„(1370) 


1200-1500 


200-500 


7171 


qq, glueball 


/o(1500) 


1505 


109 


Till 


qq, glueball 


/o(lVlO) 


1724 


137 


KK 


qq, glueball 



Table 5: Light flavour singlet 0++ mesons from the PDG. The dominant strong interaction decay is reported. 



been done in unquenched lattice QCD, and where there are possible experimental candidates. Hart 
and Teper p2| ] only found a signal for the 0++ and 2++ states from unquenched calculations. 

I summarise some pertinent experiment results for the flavour singlet 0++ mesons in table ||. 
The /o(600) (a) is particularly interesting. It was proposed in 1955, but only entered the PDG 
summary tables in 2002 [|T8|]. The parameters for the /o(600) comes from the work of Caprini et 



al. [43], while the PDG summary quotes a mass in the range 400 to 1200 MeV [29] 



Morningstar and Peardon [yTj] obtained Mo++ = 1730(50)(80) MeV for the mass of the lightest 
0++ glueball from quenched QCD. Chen et al. [0] recently found Mo++ = 1710(50)(80) MeV. 
The quark model predicts that there should only be two 0++ mesons between 1300 and 1800 MeV, 
so if the mixing between the glueball and qq operators is weak, then the 0++ glueball is hidden 
inside the /()(1370), /o(1500) and /()(1710) mesons. 

The old work by Weingarten and Lee [ p3| ] on glueball-^g mixing introduced a mixing matrix 
between the glue and qq states. The matrix elements were estimated in quenched QCD. Weingarten 
and Lee [i5 \ predicted that the /o(1710) meson was 74(10)% 0++ glueball, and hence the mixing 
between the 0++ glueball and qq states was weak. Weingarten and Lee's [ p5| ] calculations were 
critiqued in [p^]. Also their results were probably effected by the quenched artifact in the non- 
singlet 0++ correlator [19]. As discussed in section ^ now that 77, 77' mixing is now being studied 
in 2 + 1 lattice QCD calculations, then the mixing between light and strange flavour singlet scalar 
operators can in principle also be included. 

There are claims that the continuum phenomenology is more consistent with the strong mixing 
between 0++ glueball and qq states. In this case there could be a sizable contributions from glueball 
interpolating operators to the /o(600) or /o(980) mesons. 



The SESAM collaboration studied the glueball spectrum on unquenched lattices [)47[]. Mc- 
Neile and Michael studied the light 0++ spectrum with unquenched QCD [ p6| ] at a coarse lattice 
spacing and found the mass of the lightest flavour singlet 0++ meson was very light. Using 0++ 



glueball operators. Hart and Teper [ ]42| ] found that 

Mo++UNquenched = 0.85 (3)Mo++Q„enc/ierf (3-2) 

at a fixed lattice spacing of 0. 1 fm. 

Unfortunately, the existing unquenched lattice QCD calculations of the flavour singlet 0++ 
mesons don't have the range of lattice spacings where a continuum extrapolation can be attempted. 
In quenched QCD it was found that the lattice spacing dependence of the mass of the 0++ glueball 



was strong. The use of a Symanzik improved gauge action by Chen et al. [ 44 ] and, Morningstar 



and Peardon [pT]], produced a slightly smaller slope with lattice spacing of the scalar 0++ glueball 



9 



Hard Hadron spectroscopy 



Craig McNeile 



O Quenched 

<| Wilson, Wilson gauge, SESAM 

O Tad clover, Wilson gauge, UKQCD-I 

□ NP clover, Wilson gauge, UKQCD-II 
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Figure 5: Summary of unquenched results for lightest flavour singlet 0++ mesons from [|l7l|. The un- 
quenched results are from SESAM @, UKQCD-I and UKQCD-II [^. 



mass, than for calculations that used the Wilson plaquette action. This is relevant to unquenched 
calculations, because any suppression of the mass of the flavour singlet 0++ mass may be due to 
lattice spacing effects. 

The SCALAR collaboration [^], used unquenched lattice QCD, with Wilson fermions and 
the Wilson gauge action, to study the 0++ mesons. At a single lattice spacing a ~ 0.2 fm, with qq 
interpolating operators only, they obtain m-qq ~ irip . At the very least this result needs to be checked 
with a continuum extrapolation. 

In unquenched QCD, both glue and qq states will couple to singlet 0^^ mesons, so it is better 
to do a variational fit with both types of operators as basis interpolating operators. The variational 
technique analysis of the singlet 0++ mesons was done by Hart et al. [p^]. A combined fit to 
0++ glue and qq interpolating operators with two types of spatial smearing sources was done. The 
calculation used non-perturbative improved clover action at a single lattice spacing [||]. Configura- 
tions from CP-PACS [22] with the Iwasaki gauge action and tadpole improved clover action were 
also used in the analysis. A summary plot of the results is in figure ^ The data with the bursts and 
squares (with the pion masses written near them) in figure ^ shows an additional reduction of the 
mass of the 0++ state over the pure glueball operators, as used by Hart and Teper [42]. The data 
with the Iwasaki action should be less affected by lattice artifacts [22]. 

Mathur at al. [18] recently claimed to get a result for the /()(600) (a) from quenched lattice 
QCD with pion masses as low as 180 MeV. Using the interpolating operator YYsWWTsY they 
obtain niyj,(6oo) ~ 550 MeV. This interpolating operator has a disconnected contribution that they 
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computed and found to be small. The key part of this work is a three state fit {7i{p = 0)7l{p = 0) 
, /()(600), n{p = ^)7i{p = -j^) using the Bayes adaptive curve fitting algorithm [^]. They 
studied the finite volume effects to distinguish the signal for the resonance from the Tin scattering 
states [Q. 

I have a number of concerns about their result. For example, can quenched QCD really be used 
to calculate the mass of a meson with a mass of 440 MeV and width of 544 MeV p3|]? Mathur et 
al. g computed the mass of the i<:*(1430), with a width of 290 MeV to be 1 .41 ± 0. 12. 

Another concern is the use of the Bayes adaptive fitting [ ]50| ] algorithm. This is a clever way of 
doing a multi-exponential fit to a single correlator The Bayes adaptive fitting method was used by 
the Kentucky group to study the Roper resonance [|lj in quenched QCD. Other groups have not 
confirmed this (see Lasscock et al., [52] for a review), partly because the Kentucky group worked 
with much lighter pions than other groups 

Sugamuma et al. [|^] used quenched QCD to study hght 0++ states with qqqq interpolating 
operators. Hybrid boundary conditions and the Maximum Entropy Method were used to study 
correlators. No localized resonance of the form qqqq was found in the quark mass region of 
nis < niq < 2ms. 



4. Brief overview of status of unquenched QCD 



There have been individual plenary reviews for some of the recent large scale lattice QCD 



calculations []54|, p5|]. In this section I compare some results from different lattice QCD calculations. 
I used ro = 0.469 fm from HPQCD and MILC collaborations [25, to consistently determine the 
lattice spacing. I will discuss this more in section ^ This will make the comparison of some 
group's results against experiment slightly worse, but it is necessary to use a common value of vq 
to compare the results of different groups. I mostly use recently published data, rather than use the 
very latest results presented at this conference. 

I compare results at fixed lattice spacing. The results from different actions only need to agree 
in the continuum limit. In the past there was a controversy that the quenched Edinburgh plots for the 



Wilson and staggered actions did not agree in the continuum limit [57]. This was largely resolved 



by better extrapolations to the continuum limit [|58|, q9|]. All the recent large scale calculations use 
improved lattice formulations, so the corrections to the continuum limit should be small (but still 
need to be quantified of course). 



4.1 The pion decay constant 

In figure ^ I plot the pion decay constant in physical units as a function of the square of the pion 
mass in physical units. I have applied the finite size corrections from Colangelo et al. []60|]. The 



investigation of finite size effects is an active area of study at the moment (see for example []54|]). 
All the data in figure ^ either use definitions of /;;: with an automatic matching factor of 1 , or a 
renormalisation factor that is obtained via a nonperturbative method. 

I find it very encouraging that the results, in figure ^ from MILC [61] (the new results from 
the superfine run), NPLQCD [Q, HPQCD [||] and ETM collaborations [|^] cluster very closely 
together. 
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Figure 6: The pion decay constant as a function of the square of the pion mass. There is data from 
MILC NPLQCD [||], HPQCD CERN and QCDSF [M 



The data from the CERN , 65 ] and QCDSF [ 56 ] groups shghtly disagree with other groups 
at the ~ 4 MeV level. There could be several reasons for this. This could be a difference due 
to 0{a^) terms, or the way Za is treated. It would be interesting to do a similar comparison for 
heavy-light decay constants [|6^]. 

The CERN group [ 



65] presented their results as a ratio of the decay constant to a refer- 



ence pion mass. This analysis method was used by the ALPHA collaboration [59, 70 1 to look for 
chiral logs, or quenched artifacts, in the pion mass and decay constant. As previously noted by 
Jansen [58], the ratio method hides a lot of systematic errors. In figure 0, 1 use the ratio method to 



compare the data from the CERN group [|64|, |65|] with that from the ETM collaboration [|67|]. The 
agreement is much better in the ratio plot, than comparing just the decay constant. 

Now that the quality of unquenched lattice QCD calculations is improving, we may need to 
purge the subject of analysis techniques that were developed when the data was poor. For example, 
the Edinburgh (or APE) plots of the ratio of the nucleon mass to vector meson mass as a function 
of the ratio of the pion mass to vector meson mass, mix up the problems of the sensitive chiral 
extrapolation of the nucleon mass with the difficult problem of the effect of decay of the p meson. 
Usually both the vector and nucleon masses have fairly large errors, but there is no cancellation 
when the ratio is taken, so the Edinburgh plot tends to amplify the errors. 

4.2 The p meson 

In figure l8| I plot the mass of the light vector meson as a function of the square of the pion 
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Figure 7: Ratio of decay constants. Comparing the results from ETM collaboration [|7|| and the CERN 
group §51. 



mass, from lattice QCD calculations that use improved staggered, domain wall, and twisted mass 
fermions. The data in figure ^ are remarkably consistent, although the statistical errors need to be 
reduced on some results. 

In the real world the p decays into two pions, via a P-wave decay. The threshold for decay 
is 2a + where L is the side of the box. The CERN group [K4l found excited masses 



for the p channel that were consistent with lyin^ + (x)^- ^^^^ more kinematically favourable to 
study the decay of the p meson with one unit of momentum to decay to a pion at rest and a pion 



with one unit of momentum |72, 73] 



Michael and Urbach for the ETM collaboration [Q], estimated the mixing element between p 
and Tin from a three point function. This mixing produced a 5% shift in the mass of the lightest p 
using the method in ||72|]. This suggests that the mass of the p meson in figure |8| from the ETMC 
collaboration is largely unaffected by the two n decay. 

The p decay will effect the chiral extrapolation model used to extrapolate the mass of the rho 



meson. The Adelaide group have studied different regulators \ p4[ [75| ] for the effective field theory 
of p decay. This produced additional mass dependence for very light pion masses. These type 
of chiral extrapolation models for the p meson were not widely used to analyse lattice data. One 
reason for this was that the fit models showed behaviour outside simple linear or quadratic mass 
dependence in regimes where there were no lattice data. Now that many groups have unquenched 
data with pion masses in the 300 to 400 MeV regime, now is the time to revisit at this issue. Bruns 



and MeiBner [ l76[ ] have published a new chiral extrapolation formulae for the mass of the p meson. 



The derivation used a modified MS regulator and a power counting scheme. 

Mp=Ml + c,Ml + C2Ml + CiMl\n{j^) (4.1) 

The term with the C3 coefficient is due to the self energy (in the infinite volume limit). Bruns 
and MeiBner [^] recommend that the size of the c, coefficients obtained from the fits be checked 
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Figure 8: The mass of the light vector meson as a function of the square of the pion mass. I include prelim- 
inary data from ETMC and pubhshed data from RBC-UKQCD and the MILC collaboration 




Figure 9: Quark mass dependence of the mass of the p 



against other estimates (eg. from (Opn coupling). Figure g from the paper by Bruns and MeiBner [ ]76[ ] 



shows the effect of the extrapolation model in equation |4JJ on older CP-PACS data. See | |77| ] for 
an another analysis of the p decay on the mass of p meson. 



5. Determining the lattice spacing 

To convert the results from unquenched lattice calculations from lattice units into physical 
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units, requires that the lattice spacing is determined. In principle this requires that a single physical 
number be sacrificed to determine the lattice spacing. There are problems with many of the obvious 
choices of quantities to use. For example the mass of the nucleon is very sensitive to the chiral 
extrapolations of the mass to physical quark masses. Now that most groups are using some kind of 
non-perturbative renormalisation in the calculation of the fj^ this is in principle a good quantity to 
use to determine the lattice spacing. Unfortunately, the determination of fji also requires a careful 
chiral extrapolation. Also as Marciano [ f/Sl ] pointed out the ratio ^ can be used to extract the CKM 
matrix element Vug. 

Quantities that are not good choices to determine the lattice spacing are hadrons that decay via 



the strong interaction [79]. Given the discussion about quark mass dependence of the mass of the 



hght vector meson in section 4.2, it is clear that using the p mass or K*' mass to set the scale is not 
a good idea, because we don't have full theoretical control of the chiral extrapolation. The problem 
of using the p mass to set the scale has been known for a long time, for example it was discussed 
by DeGrand at the lattice 90 conference [pO[], although it is still being used J/I]]. 

However, I don't think we should give up studying the properties p or K* mesons using lattice 
QCD. Obtaining an accurate value of the p mass from lattice QCD is prerequisite to calculations 
of hadrons such as the Roper resonance or scalar mesons. There are also a number of important 
semi-leptonic weak decays, such as B p + y, B + y that involve a p meson in the final 
state. 



One popular way to determine the lattice spacing is to use a quantity called ro 081|], determined 
from the heavy quark potential, via 

F{ro)rl = 1.65 (5.1) 

where F is the force of the heavy quark potential. The value of ro in lattice units can be obtained 
very accurately from lattice calculations. When Sommer [^] introduced ro as he assigned ro = 
0.49 fm, from continuum heavy quark potentials, with an estimated uncertainty of 10%. The value 
of ro has largely replaced the use of the string tension to determine the lattice spacing. 

The MILC/HPQCD collaborations have determined ro = 0.469(7) fm from the 2S-1S mass 
splitting in T from unquenched lattice QCD calculations with 2+ 1 flavours of improved staggered 
quarks [p6|]. A consistent lattice spacing is determined from other quantities, such as other splittings 



in the T system, fj^, and mass splittings in charmonium [J79[]. This value of ro is crucial to the 
phenomenology produced by the improved staggered program. For example, if I naively change 
from using ro = 0.469 fm to ro = 0.49 fm, with HPQCD's recent result for fo, = 241(3) MeV, then 
their "new" number is 231(3) MeV [fl. This is a shift of 3 a. As the HPQCD collaboration push 



to ever higher precision the current error on ro needs to be reduced [|63|]. 

Other collaborations, for example RBC-UKQCD [^T]], use ro = 0.49 fm to set the lattice spac- 
ing, based on consistency with p mass and the K*/K mass ratio. As reviewed by Boyle []55|], the 
RBC-UKQCD collaboration are now obtaining lower values of ro than 0.49 fm, when they set the 
lattice spacing using hadrons that are stable under the strong interactions, such as the Q.. 

The value of ro can be determined from lattice calculations. In table ^ I collect some estimates 
of the value of ro from lattice QCD calculations. From the summary of the results in table ^ we see 
that the recent unquenched calculations are starting to report values of ro between 0.44 to 0.48 fm. 
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Group 


Uf 


Method 


ro fm 


Sommer [^] 
Morningstar and Peardon [p2|] 





quark model 
quenched summary 


0.49^^' 
0.48(2) 


UKQCD [§] 
JLQCD Q 
JLQCD Q 
QCDSF [|4|] 


2 
2 

2 


K/K* 
ntp 

summary nucleon masses 


0.55 
0.497 (-9)(13) 
0.5702(75)(50) 
0.47(3) 


QCDSF[77| 
ETMC [67] 
HPQCD/MILC/FNAL ^ ||] 


2 
2 

2+1 


k 

gA 

fn 

Upsilon & ratio plot 


0.45(1) 
0.454(7) 
0.469(7) 



Table 6: Summary of lattice determinations of tq. 



I am being slightly hypocritical here, because for clover action with niq > rus/l we used to 



argue that by using ro ~ 0.49 fm some systematics might cancel [4C]. 



One problem with using rg to set the scale is that rg needs a chiral extrapolation to this massless 



limit. See Sommer et al. [|850, Aoki [ ]57| ] and Bhattacharya et al. [ |86p for discussions about the mass 
dependence of ro. Some part of the mass dependence of ro may be physical, and some part may be 
a lattice artifact. In figure |lO|l plot ro, with the intercept fixed at 0.467 fm, for a number of different 
lattice formalisms, against the square of the pion mass in physical units. 

The clover action is only an onshell improved action. In the ALPHA formalism for the im- 
provement of the clover action, the leading lattice artifact part of the linear mass dependence of ro 
is related to the bg improvement coefficient. 

gl = {l+bgamci)gl (5.2) 

Fully 0{a) improved actions, such as domain wall or twisted mass fermions should have no 
0(a niq) lattice artifacts, and reduced dependence on the quark mass on ro. I note that RBC- 
UKQCD [j8^, ^ linearly extrapolate ro with quark mass, but the ETM collaboration [ p7| ] extrapo- 
late ro quadratically in quark mass to the chiral limit. 

My naive understanding of the phrase "decoupling of continuum and chiral limits" about 
fermion operators that are approximate solutions to the Ginsparg-Wilson relation suggested that 



mass dependence of ro for overlap/DWF actions should be small. The data in figure [10| doesn't 
appear to show reduced mass dependence for the domain wall data. The domain wall data is at a 
coarser lattice spacing than for the other actions, so I may be making an unfair comparison. 



The MILC collaboration [gjl have chosen to work at fixed lattice spacing rather than fixed j8, 
so I have not included their results in figure |l^. As the quark mass is changed, jS is changed to 
keep ro/a constant. This was also the strategy of the improved clover calculations performed by 



the UKQCD collaboration [§, g9p. 

In this section I am not advocating that ro should always be used to set the lattice spacing in 
every calculation. However, I do hope that eventually there will be consensus between different 
lattice QCD calculations on a final value for ro in physical units. It is certainly a crucial cross-check 
on the results from the improved staggered program. 
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Figure 10: The value of ro as a function of the square of the pion mass, using data from RBC-UKQCD ||7T[], 
QCDSF-UKQCD [M, and ETMC 



6. Conclusions 



There is still no consensus as to whether qq operators in lattice QCD calculations are cou- 
pUng to the fl'o(980) meson. To clear up the many questions about the spectrum of the 0++ scalar 
mesons, unquenched lattice QCD calculations with tetraquark interpolating operators are required. 
There is "some" evidence that the flavour singlet 0++ interpolating operators, in unquenched lat- 
tice QCD calculations, are coupling to states around or below a 1 GeV [4^9]. Although a continuum 
extrapolation is required for definite results and the open decay channel issue needs to be studied. 

After all the successful work on algorithms for reducing the mass of the sea quarks in lattice 
QCD calculations, we are now starting to study light mesons with open S-wave decay channels. I 
presented some evidence from the ETM and MILC collaborations, that the bi meson has an open 
strong decay channel. This type of issue will be of increasing importance for lattice studies of 
particles that decay via the strong force. For example, the MILC collaboration [ ]90| ] claimed to see 
problems with the light exotic 1^+ meson, because of an open decay channel. Eventually, the issue 
of dealing with resonances in lattice QCD will be dealt with by Luscher's technique [|9T|]. Until 
then pragmatic approaches to studying strong decays on the lattice are still important. This year 
Luscher's technique for resonances was applied to the p meson for the first time, by the CP-PACS 
collaboration [|73[]. 



In the past few years the MILC collaboration [|24|, ^] have been doing the unquenched calcu- 
lations with the lightest pions, largest physical volumes, biggest range of lattice spacing. This year 
mostly due to algorithm improvements, and partly due to bigger machines other collaborations now 
have lattice QCD results that are comparable to the quality of those from the MILC collaboration. 
It is important to do cross-checks on the results from different lattice QCD formalisms, as part of 
the quest for errors at the percent level. Many of the older ways to compare results from different 
lattice QCD calculations may no longer be appropriate to this high precision era of lattice QCD. 
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